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Abstract 

The results of evaluation of mixed QEDxQCD corrections to R{s) in e + e~ 
annihilation and R T in hadronic decay of the r lepton to 0(a n a™), m + n < 3, 
are presented. The net effect on a s (Mz) from the Z decay is only about 0.1% 
and in the r decay case the net effect is negligible, as expected. 
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Inclusive processes of the hadronic decays of the Z boson and the r lepton are in the 
center of experimental and theoretical |||4[] considerations. These processes are best 
suitable for, e.g, high and low energy determinations of the strong coupling a s (s), thus 
providing one of the crucial tests of the Standard Model (SM). High experimental and 
theoretical precisions are important to to spot any deviations from the SM and thus to 
find a room for a possible new physics. The precise a s (s) is also crucial for checking of 
fundamentals of GUT B. The up to date results for a s extracted at different energies, using 
various observables are in satisfactory agreement with the SM ||. However, the problem 
of unification of couplings at the GUT scale remains [|5|. There has been much attention 
drown to the so called Rb — R c problem ||. In the above examples, high experimental 
and theoretical precisions are necessary. A significant progress has been made in the last 
decade in theoretical evaluation of the above quantities, using perturbation theory methods 
To resolve the observed discrepancy between the experiment and the SM predictions, 
numerous calculations of contributions beyond the SM has been done recently (see, e.g., ||). 

In the present work we evaluate corrections due to multiple photon exchange between 
quarks and between leptons taking place along the gluon exchange between "final state" 
quarks. The calculation has been done to four-loop level corresponding to 0(a n a™) (m+n < 
3) corrections. The objective of this work is to check explicitely that there are no anomalously 
large perturbative coefficients in front of a small a n a™ expansion parameter, that could make 
the correction significant. Such a possibility is not excluded (see, e.g., 0). Indeed, individual 
diagram (or even gauge invariant sets of diagrams) contributions are often large and only 
thanks to some delicate cancellations, they add up to small numbers. 

We found that all calculated new corrections are indeed small and, as expected, they can 
be neglected at the relevant level of precision of 0(1%). 

The method of calculation is similar to the one used in earlier works riXTXI and described 



in ref. ||11|| . Here we discuss issues specific mainly for the present calculation. Known 
electroweak and QCD corrections to the Z decay rates are summarised in refs. || and [0, 
corresp ondingly. 

The decay rate of the Z boson is usually calculated as the imaginary part of the correlation 
function of two neutral weak currents of quarks coupled to Z boson. The decay rate has 
vector and axial- vector parts, which generally are different. In the limit of massless quarks - 
the approximation used in this work, the vector and axial parts are identical up to the known 
factors. Therefore, here we calculate only a nonsinglet vector part, or the quantity R(s). 
The total hadronic decay rate V z can then be obtained straightforwardly, by multiplying the 
result on the sum of the squares of vector and axial couplings and by adding other known 
corrections from refs. PJTJ. 

The four-loop R(s), including the mixed QCDxQED corrections of 0(a n a™) (m + n < 3) 
has the following form 
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The powers of electromagnetic coupling a and the strong coupling a s in the superscripts 
of Ri indicate correspondingly the number of photon and the number of gluon exchanges 
involved in the corresponding multiloop Feynman graphs. The couplings are renormalized 
according to the MS prescription |T2] 
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where the couplings run according to the renormalization group equations 



a(3 QED (a,a s ) = /r 
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and the two-loop renormalization group (3 functions now include 0(aa s ) perturbative terms 
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The known values for one- and two- loop /3 function coefficients can be found, e.g., in ref. 
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The 0(aa s ) coefficients can be found from the analysis of two- loop graphs and the known 
two-loop results 
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p{a 3 a) = —~T Yl Q% (3 { r s) = -\c F N F E Q). (7) 

f=u,d,... f=u,d,... 

The eigenvalues of the Casimir operators for the adjoint (Na = 8) and the fundamental 
(Np = 3) representations of SU C (3) gauge group are Ca = 3 and Cp = 4/3. The Dynkin 
index T for the fundamental representation is usually chosen T = 1/2. Ni is the number 
of leptons appearing in closed fermion loops not attached to the gluon propagators. Qf 
denotes electric charge of the quark of flavor / in the units of the electron charge. 

In order to calculate the unknown perturbative terms in Eq. (|I|), we have rederived Eq. 
(6.5) of ref. [[□]] - the four-loop expression for R(s) in terms of perturbative coefficients 



of the correlation function, its renormalization constant and the (3 function, taking into 
account 0(a m a™) terms. We have also found the renormalization group constraints, similar 



to the ones in Eqs. (2.29), (2.30) and (2.32) of ref. [11]. These relations are very helpful in 
testing the intermediate results. Contributions from two-, three- and four-loop graphs were 
found using the graph- by-graph results of the work [|10j. This requred a careful counting 
of symmetric and statistical factors for each of the contributing 98 four-loop, 14 three-loop 
and 2 two-loop Feynman graphs and calculation of the gauge group weights. 

We obtain the following analytical results for the QED and QEDxQCD terms in Eq. ([!]) 
for the standard U(l) and SU C (3) gauge groups within the MS framework 

R{ a) = 9 jY,Q% (8) 
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where summations run over participating quark flavors, normally - u, d, s, c, b. We keep 
summation indices different in order to identify its source topology. For instance, the sum- 
mation over / is due to the outer quark loops - the "final state" quarks. The summations 
over i and j correspond to the inner quark loops where virtual quarks and leptons can ap- 
pear. Note that, of course, leptons can appear only in the case when there is no gluon line 
attached to the fermion loop. The last terms in Eqs. (pi]) and ([H]) are due to the specific 
four-loop topology where current operators are inserted in the separate fermion loops. The 
two-, three- and four-loop QCD terms are known and can be found in the original work |jTU . 
Corrections due to non-vanishing quark masses to these terms are also known and can be 
found, for instance, in the review paper |7|. For simplicity, in the present work we neglect 
all mass corrections which can trivially be added to our results. 

For five massless quarks and infinitely heavy top quark we obtain the following numerical 
result 
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One can see that the calculated corrections are very small. The largest corrections beyond 
the leading QED term - ~ a 2 , ~ a s a and ~ a 2 a have a same sign and they add up to the 
total effect of about 0.1 MeV on the Z width and only 0.1% on a s . The other calculated 
corrections are clearly negligible for any phenomenological applications. The calculated 
terms of 0(a 2 ) and 0(aa s ) confirm known results |y| . 

Let us use the calculated corrections on the Z width and obtain the similar terms for the 
quantity 
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the perturbative part as an integral (see, e.g., [|Tl|] for details) 



in the r lepton hadronic decay process. Following the well known procedure |I4| , we calculate 
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where IT(s) is the transverse part of the correlation function of weak currents of quarks 
coupled to W boson. The countour C is the circle of radius \s\ — M 2 . In the integrand we 
substitute the a s (s) and a(s) by their expansion in terms of a s (M T ) and a(M T ) using the 
equations 
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Then we evaluate contour integral and express R T in terms of the perturbative coefficients Ri 
in Eq. (Q), making obvious substitution X)/ Q 2 — * \Vud\ 2 + \Vus\ 2 and omitting terms coming 
from the four-loop topologies, where current insertions are placed in the separate fermion 
loops. Thus for the perturbative part of the R T we obtain the following MS analytical result 



RT = Z(\V ud \ 2 + \V US | 2 ){l + 



7T 



7T 



a(M r )y_ (Q 3 ) 



| a s (M T ) ^( Qs) + a(M r ) r ( aas) + 



7T 



7T 



7T 



7T 



T .(« 2 « S ) 



+ 



a s (M T ) 



7T 



Jo?) + 



a(M T ) 



+ 



7T 



7T 



r 3 



(19) 



where for the standard U(l) and SU C (3) gauge groups within the MS framework we obtain 
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The summation index / runs over w, d, s "final state" quark flavors appearing in the outer 
fermionic loops and the summation over j corresponds inner (virtual) fermionic loops. We 
take j = u, d, s. For the number of virtual leptons in the inner fermionic loops, we will take 
Ni = 2 (e, /i). In fact, the r-lepton can also appear in the virtual loops. However, within 
our approximation, masses of all quarks and leptons, including virtual ones are neglected. 
On the other hand, the r mass cannot be neglected at this scale. Therefore, the effect of 
the virtual r lepton is left out of our consideration. This effect can be evaluated via exact 
evaluation of Feynman graphs with massive loops. The corresponding numerical effect is 
expected to be small. The two-, three- and four-loop QCD terms are known and can be 
found in [jig] . We obtain the following numerical expression 
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As one can see, the higher order QED and QEDxQCD corrections are very small. Even 
the largest ones of 0(aa s ) and 0(aa 2 s ) have similar size but opposite signs and they cancel 
each other along the 0(a 2 ) correction. Thus, what remains is the leading 0(a) correction 
that has some numerical relevance. All other corrections are clearly negligible. 

Summarising, we note that the calculated higher order QED and mixed QEDxQCD 
corrections to the hadronic decay rates of the Z boson and the r lepton can be safely 
neglected in present phenomenological applications. 
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